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(54) Metfiod, apparatus and computer program product for forming data to be analyzed by finite 
element method and calculation method based on finite element method 

(57) A user uses part data forming means 101 and 
module data forming means 102 to input the actual 
dimension, physical constants and mesh dividing 
number for fundamental shapes which are registered in 
advance, thereby forming parts, and then indicates the 
relative position between the parts to form the entire 
shape of an assemt)ly of plural parts without paying 
attention to coincidence or non-coincidence of nodal 
points. Data converting means 104 divides the shape of 
each part thus assembled according to the indicated 
mesh divisional number to generate element data and 
nodal point data. Further, it generates a constraint 
equation for connecting nodal points which are non- 
connected between neighboring parts, and forms an 
analysis model 401. A finite element method analyzer 
1 05 uses approximate calculation means 1 06 to approx- 
imate a non-connected nodal point displacement from a 
nodal point displacement of neighboring structural ele- 
ments on the basis of the constraint equation. 
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Description 

BACKGROUND OF THE INVENTION 

Field of the Invention: 

The present invention relates to method, apparatus, 
and computer program product for forming data to be 
analyzed by finite element method (FEM) by which dis- 
placement, temperature, and the like are analyzed and 
a calculation method based on finite element method. 

Description of the Prior Art: 

According to the finite element method, the shape 
of an object which is an analysis target is divided into 
small units, a solution is sectionally approximated by a 
function having a finite value for every units to calculate 
the solution of the entire shape. Hereinafter, the division 
of the entire shape Into small units is referred to as 
"mesh division'', each unit thus obtained to as "ele- 
ment", and an apex of an element to as "nodal point". 
The shape to be analyzed is defined by an assembly of 
element data and nodal point data, and these data are 
refen'ed to as "structural data". 

In order to the form structural data, it is necessary 
to form the overall shape to be analyzed. A method of 
inputting the coordinates of all the apexes of the shape 
has been hitherto widely used to form the overall shape. 
However, this method has the disadvantage that it 
needs more time as an analysis target becomes three- 
dimensional and more complicated. 

Therefore, JPA-6-4630 disclosed a method of com- 
bining fundamental shapes and local shapes which are 
registered in advance, thereby forming the overall 
shape. According to this method, the basic shapes and 
the local shapes are selected to input the actual dimen- 
sions, and the respective shapes of the basic shapes 
and local shapes are combined with one another, 
whereby the overall shape can be readily formed. 

In order to form the structure data, a mesh dividing 
work is further needed. Ttiere have been hitherto known 
various techniques for supporting the mesh dividing 
work. For example, according to the method disclosed 
in JPA-6-4630, the overall shape of the assembled ele- 
ments is automatrcally divkled into hexahedral regions, 
such divisional numbers in the longitudinal, lateral and 
height directions that are equal between neighboring 
hexahedral regions are received from a user for every 
hexahedral regions, and the mesh division is performed 
on the basis of the divisional numbers. Further, accord- 
ing to the method disclosed in JPA-3-29057. the mesh 
division is automatically performed by a three-dimen- 
sional solid mesh generator, and each nodal point is 
moved by using a three-dimensional isoparametric 
smoothing method so that the shape of each element is 
well regulated. Still further, according to the method dis- 
closed in JPA-5-108694, a similar shape model which is 



similar to a shape model to be analyzed and formed as 
an assembly of square blocks is generated, and it is pro- 
jected onto the shape model (analysis target) to perform 
the mesh division. Any technique as explained above is 

5 based on the idea that the nodal points of the respective 
elements constituting the analysis target shape are 
coincident between neighboring elements. 

As described above, it has been a dominant idea in 
the conventional techniques that the nodal points must 

10 be coincident between the neighboring elements. 
Therefore, when some fundamental shapes are com- 
bined to form the overall shape as disclosed in JPA6- 
4630. much time is needed for the mesh dividing work 
thereof. This is because a user must indicate the mesh 

IS divisional number so that the nodal points are coinci- 
dent between neighboring areas. 

gUMMARYQFTHE INVENTION 

20 Therefore, an object of the present invention is to 
provide a method for forming data to be analyzed by a 
finite element method, which can perform the work of 
forming the shape to be analyzed and of mesh division 
in short time. 

25 Another object of the present invention is to provide 
a method for executing calculation based on a finite ele- 
ment method which can obtain a solution even when the 
nodal points are not coincident between elements. 
According to one aspect of the present invention, 

30 there is provided a method for forming structural data to 
be analyzed by a finite element method which com- 
prises steps of: providing structural elements having 
various shapes; inputting dimension for each of the 
structural elements; applying mesh generation to each 

35 of the structural elements, and piling up the structural 
elements irrespective of coincidence or non-coinci- 
dence of nodal points. 

According to another aspect of the present inven- 
tion, there is provided a calculation method based on 

40 finite element method for calculating nodal point dis- 
placement or nodal point temperature In case that struc- 
tural data to be analyzed by a finite element method 
contain non-connected nodal points, which comprises a 
step of calculating the displacements or temperatures of 

46 non-connected nodal points of a structural element by 
using the displacements or temperatures of nodal 
points of the neighboring structural elements on the 
basis of constraint equations. 

According to still another aspect of the present 

50 invention, there is provided an apparatus for generating 
input data to be analyzed by a finite element method, 
which comprises: a storage device which stores funda- 
mental shapes; part data forming means for setting the 
actual dimension, the physical constants and the mesh 

55 dMding number for the fundamental shapes in accord- 
ance with a user's instruction, thereby forming part data 
for the fundamental shapes; module data forming 
means for forming the data of a module on the basis of 
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a pair of parts indicated by the user and the relative 
position between the parts in the pair, the module being 
constructed by a combination of the paired parts; analy- 
sis condition data forming means for forming analysis 
condition data according to the analysis condition Indi- 
cated by the user; and data converting means for input- 
ting model data comprising the part data, the module 
data and the analysis condition data, and forming an 
analysis model comprising structural data, converted 
analysis condition data and constraint equations, 
wherein the structural data comprises elements data 
and nodal points data which are obtained by building up 
the parts in accordance with the relative positions and 
then dividing the shapes of the parts by the mesh divi- 
sion number, each components of the converted analy- 
sis condition data corresponds to each nodal points, 
and the constraint equations connect non-connected 
nodal points between the parts. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a flowchart showing an example of the 
process flow of a data forming method fa a finite 
element method according to the present invention; 
Fig. 2 is a diagram showing an exanple of struc- 
tural data obtained by piling up a cylindrical (colum- 
nar) structural element on a rectangular 
paralletopiped structural element; 
Fig. 3 is a diagram showing an example of a con- 
nection status of a nodal point non-coincidence ele- 
ment of a structural element connection portion; 
Fig. 4 is a diagram showing a rectangular parallele- 
piped structural element; 

Fig. 5 is a diagram showing a triangular prism struc- 
tural element; 

Fig. 6 is a diagram showing a cylindrical structural 
element; 

Fig. 7 is a spherical structural element; 

Fig. 8 is a diagram showing a trapezoidal hexahe- 

dral structure element; 

Fig. 9 is a diagram showing a notched pillar struc- 
tural element: 

Fig. 10 is a square cylindrical structure element; 
Fig. 11 is a table showing input items which are 
needed to define a figure; 
Fig. 12 is a table showing a comparison example of 
the number of numerical values to be input when a 
structural element is defined, and the number of 
numerical values which is needed to define a figure 
on the basis of a conventional coordinate input 
method; 

Fig. 13 is a diagram showing a structural element 
assembly example; 

Fig. 14 is a table showing an example of a input dis- 
play frame when a silicon rectangular parallelopi- 
ped structural element is assembled with an 
aluminum rectangular parallelepiped structural ele- 
ment set as a base; 



Rg. 15 is a diagram showing a constraint condition 
set; 

Fig. 16 is a table showing a comparison example of 
the minimum value and the maximum value of the 
5 maximum main stress of models when nodal points 
are continuous and when nodal points are discon- 
tinuous; 

Rg. 17 is a block diagram showing an embodiment 

of a data forming apparatus for a finite element 
10 method according to the present invention; 

Rg. 18 is a flowchart showing an example of the 

processing of a part data forming means; 

Rg. 19 is a diagram showing part data; 

Rg. 20 is a flowchart showing an example of the 
75 processing of a module data forming means; 

Fig. 21 is a diagram showing module data; 

Rg. 22 is a diagram showing analysis condition 

data; 

Rg. 23 is a flowchart showing an example of the 
20 processing of a data conversion means; 

Fig. 24 is a diagram showing an analysis model; 
Fig. 25 is a diagram showing a connection status of 
a nodal point non-coincidence element of a part 
connection part; and 
25 Rg. 26 is a diagram showing ^,ti,^ normal coordi- 
nates. 

DETAILED DESCRIPTION OF THE PREFERRED 
EMBODIMENTS 

30 

[embodiment] 

Preferred embodiments of a data forming method 
and a calculation method for a finite element method 
35 according to the present invention will be explained 
hereunder with reference to the accompanying draw- 
ings. 

Rg. 1 is a flowchart showing the process flow of a 
data forming method for a finite element method 
40 according to the present invention. This process is per- 
formed on a conputer system which includes a compu- 
ter such as a work station or the like, a storage device 
such as a magnetic disk or the like, an input device such 
as a keyboard, a mouse or the like and a display device 
45 such as a CRT display or the like. Some structural ele- 
ments having simple shapes are registered in the stor- 
age device, and users form structural data using these 
structural elements. 

Rrst, a user inputs dimension data and data of 
so physical constants for all the structural elements which 
are usable for a final target shape (SI, S2). and then 
performs a mesh dividing work before assembling these 
structural elements (S3). Each of the structure elements 
divided into meshes can be stored as structural element 
55 DX in the unit of structure element (S4). The order of the 
Input of the dimension data and the physical constants 
and the mesh dividing work is not limited to the above 
order. In this embodiment, since the mesh division is 
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carried out on a structural element having a simple 
shape, automatic mesh generation used In a conven- 
tional finite element method is not necessary and the 
mesh division of each structural element Is performed 
by inputting a division number for the sides of the struc- 5 
tural element or the like. For example, in the case of a 
rectangular parallelepiped, a divisional number is input 
for both sides in longitudinal and lateral directions of a 
base and for a side in the height direction to perform the 
mesh division. In the case of a cylinder, a divisional 10 
number is input for each of the circumferentia). radial 
and height directions to form the mesh division. 

When the input of the dimensions and the physical 
constants for necessary structural elements and the 
mesh division are completed (YES: S5). these structural is 
elements are assembled to construct structural data. 
First, a structural element serving as a base is selected, 
and the relative position relative distance of any apex of 
a next structural element with respect to any apex of the 
base structural element set as an origin and rotating 20 
angle of the next structural element are irput to pile up 
structural elements (S6). The pile-up result can be sep- 
arately stored as assembled structural data DY (S7). 
The structural elements are piled up as explained above 
to thereby construct structural data to be analyzed 2s 
(analysis target). At this time, it is sufficient to select any 
apex of any piled-up structural element as the origin and 
then input any apex coordinate of a structural element 
to be piled up with respect to the origin. Therefore, it is 
unnecessary to take coincidence or non-coincidence of so 
nodal points between the structural elements into con- 
sideration. Further, in this structural data forming 
method, not only the structural element data DX stored 
are usable, but also the structural data DY which have 
been asserrfbled are usable like the structural element 35 
dataDX. 

After the structure data to be analyzed has been 
constructed (YES: S8), the data are con-ected if neces- 
sary (S9). In the data correction, not only the pile-up 
position is changed, but also the dimension, the physi- 40 
cal constants and the mesh division may be changed for 
every structural elements. 

Fig. 2 is a diagram showing an example of struc- 
tural data when cylindrical structural element 1 is piled 
up on rectangular parallelepiped structural element 2 in 45 
the process as explained above. As shown in Rg. 2, the 
nodal points of the elements on the upper surface of 
rectangular parallelepiped structural element 2 are not 
coincident with the nodal points of the elements on the 
lower surface of cylindrical structural element 1 . so 

Next, a method of calculating nodal point displace- 
ment or nodal point tenperature when the structural 
data contain non-connected nodal points will be 
explained. 

Fig. 3 shows a connection status of nodal point ss 
non-coincidence elements at a structural element con- 
nection portion. In Fig. 3, reference numeral 4 repre- 
sents an element of a structural element, and it is 



designed in a rectangular parallelepiped shape and has 
nodal points LULLaaaQ. Reference numeral 3 rep- 
resents an element of another structural element which 
is piled up en element 4. and m represents a nodal point 
located on the upper surface of element 4. Fig. 3 shows 
an example of the connection status that the rectangu- 
lar parallelepiped elements are connected to each 
other, and no limitation is imposed upon the shape and 
dimension of the elements. 

A displacement vector at the non-connected nodal 
point OQi of element 3 in Fig. 3 is approximately calcu- 
lated by the following constraint equation which uses 
the displacements of the four nodal points 1. i, k. ! en the 
same plane of element 4 and interpdative coefficients: 



mspnterpolative function] 



(1) 



In equation (1), the displacement of non-connected 
nodal point m is interpelatively calculated from the dis- 
placements of the nodal points at the four points en the 
same plane of neighboring element 4, and no restriction 
is imposed upon the number and position of nodal 
points for use. For example, the displacement of non- 
connected nodal point m may be interpelatively calcu- 
lated using the displacements of all the nodal points L 1. 
Is, L 0. a a Q of neighboring element 4 in Fig.3. 

When the structural data contain non-connected 
nodal points, the displacements off all the non-con- 
nected nodal points are approximately calculated by the 
constraint equation using the interpolative function. Fur- 
ther, when another structural element is piled up on the 
structural element comprising an assembly of this kind 
of elements, the same calculation as explained above is 
successively performed to thereby obtain all the resolu- 
tions. 

In the themial conduction analysis using nodal 
point temperature in place of the nodal point displace- 
ment, by using the temperature values T|, Tj, T|^, T| of 
nodal points L i, k, 1 of element 4 in place of the displace- 
ment vectors of element 4 of equation (1). temperature 
Tfn of the non-connected nodal point m of element 3 can 
be also interpelatively calculated in the same manner as 
explained above. In the interpolative calculation of the 
non-connected nodal point temperature, no restriction is 
imposed upon the number and position of nodal points 
for use. 

No restriction is imposed upon the interpolative 
function for use. however, an excellent approximate 
value can be obtained when an isoparametric shape 
function is used. 
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[example] 

Next, an example of the data forming method and 
the calculation method for the finite element method 
according to the present Invention will be explained in 
detail. 

No restriction is imposed upon the type and number 
of the structural elements which are registered in 
advance. However, in this example, it is assumed that 
seven structural elements shown In Figs. 4 to 10 are 
registered in advance and an application analysis tool is 
formed for evaluation. Fig. 11 tabulates input items 
needed to define figures shown in Figs. 4 to 10. 

In the case of a triangular prism (Fig. 5). a trapezoi- 
dal hexahedron (Fig. 8). a notched pillar (Fig. 9). the 
divisional input number is set to as small a value as pos- 
sible in order to simplify the divisional work. That is, for 
the triangular prism, (divisional number of side C) = 
(divisional number of side B), for the trapezoidal hexa- 
hedron, (divisional number of side B) » (divisional 
number of side A) and (divisional number of side C) = 
(divisional number of side C), and for notched pillar, 
(divisional number of side A) + (divisional number of 
side B) = (divisional number of notched arc) and (divi- 
sional number of side d') = (divisional number of side d) 

Fig. 12 tabulates a comparison result of the number 
of numerical values to be input in order to define these 
structural elements between this example and the con- 
ventional coordinate inputting method. However, in the 
conventional method, when the mesh divisional number 
is required to be partially input, the input number more- 
over increases by the number required to define the 
mesh divisional number. It is assumed that for a cylindri- 
cal shape, (number of apexes) = (divisional number of 
circumference) x 2. and for a spherical shape, (number 
of apexes) ^ (divisional number of circumference) x 
(divisional number of circumference ■ 2) + 2. and (divi- 
sional number of circumference) = 36 in Fig. 12. In the 
conventional method, it takes much time to calculate the 
coordinates of a structural element. On the other hand, 
in this example, the input of figures' data can be per- 
formed in short time because the dimensions of the fig- 
ures are directly input. Further, in the conventional 
method, as the structure is more complicated, it is more 
difficult to determine the coordinates. On the other 
hand, in this example, the structural element is defined 
by only the dimension thereof, and the position thereof 
is defined by setting the coordinate of any one apex 
thereof with respect to an origin which is set to any apex 
of any other structural element. Therefore, even if the 
structure is complicated, the input is still simple. 

Fig. 13 shows an example of an assembly of struc- 
ture elements. Aluminum square cylindrical structure 
element 6 and silicon rectangular structural element 8 
are disposed on aluminum rectangular parallelepiped 
structural element 5 set as a base, and then copper 
square cylindrical structural element 7 is disposed on 
aluminum square structural element 6. In this case, 



each of alt the structural elements other than silicon rec- 
tangular parallelepiped structural element 8 has an 
apex which is coincident with that of the structural ele- 
ment serving as the base, and thus the calculation of 

5 the relative position thereof is unnecessary because 
when a new structural element is added, any apex of 
any already built up structural element can be set as 
base point. For silicon rectangular parallelepiped struc- 
tural element 8, it is sufficient only to take the coordinate 

10 of one point (apex 5" in the figure) into consideration. In 
addition, tiie origin serving as tiie base is set to an apex 
(apex 1' in the figure) on the upper surface of aluminum 
rectangular parallelopiped structural element 5, and the 
relative distances of apex 5" with respect to apex V 

75 along x-direction and y-direction are merely input. 
Therefore, tiie Input time can be greatly shortened as 
compared with the conventional method in which all tiie 
apex coordinates of the structural body to be analyzed 
are input. 

20 Rg. 1 4 shows an example of an input display frame 
when silicon rectangular parallelopiped structural ele- 
ment 8 is disposed on aluminum rectangular parallelopi- 
ped structural element 5 set as a base. T3. PAT 
represents an identifier for silicon rectangular parallelo- 

25 piped structural element 8. In this case, silicon rectan- 
gular parallelopiped structural element 8 alone is 
disposed on aluminum rectangular parallelopiped struc- 
tural element 5, and thus both of "member name of 
assembly target" and "part name of origin of member" 

30 are set to T3. PAT 

However, when an assembly of plural sti'uctural ele- 
ments which have been assembled is piled up on 
another structural element, an identifier for the assem- 
bly is displayed in "member part of assembly target", 

35 and an identifier for a structural element which is in the 
assembly and serving as a base of these elements is 
displayed in "part name of origin of men*ers". T1.PAT 
represents an identifier for aluminum rectangular paral- 
lelopiped structural element 5. and displayed in "part 

40 name of starting point of assembly". In Fig. 14, tiie rela- 
tive distance (1 5, 1 3, 0) is input while apex 5" of the bot- 
tom surface of silicon rectangular parallelopiped 
structural element 8 is set to "origin number" and apex 
V on the upper surface of aluminum rectangular paral- 

45 lelopiped structural element 5 is set as the starting 
point. The rotational angle is set to zero. Since an initial 
value is equal to zero, it is practically sufficient to input 
the relative distance between the starting point and the 
origin along only X-axis and Y-axis. 

so The thermal stress when the assembled structural 
body is varied in temperature from 25''C to 760*'C is cal- 
culated by the finite element method. As shown in Rg. 
1 5, the constraint condition is set so that one apex of tiie 
bottom surface of aluminum rectangular parallelopiped 

55 structural element 5 is constrained in the x-. y-, and z- 
axis directions, the diagonal apex is consti'ained in ttie 
z-axis direction, and one of the residual apexes is con- 
strained in ttie X-, and z- axis directions while the ottier 
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is constrained in the y-, and z- axis directions. A con- 
straint equation based on an isoparametric function is 
used for the irtterpolative calculation of the non-con- 
nected nodal point displacement. In a preproduction 
tool, the meshing work is performed before assembly, 5 
however, a nodal-point-coincidence model is also pre- 
pared for a calculation precision test for nodal-polnt- 
non-cotncidence. Fig. 16 shows the minimum value and 
the maximum value of the maximum main stress for 
nodal-points-continuous model and nodal-point-discon- 10 
tinuous model. The difference of the maximum main 
stress between the nodal-point-continuous model and 
the nodal-point-discontinuous model is within a tolera- 
ble range from the viewpoint of simply comparing the 
difference in the structure and material. If a model in 75 
which the divisional number is set at a large value is 
formed, the distance between nodal points is reduced, 
and also the precision of the isoparametric interpolation 
is improved, so that the difference between the nodal- 
point-continuous mode and the nodal-point-discontinu- so 
ous model is further reduced. 

[embodiment] 

Next, an embodiment of a data forming apparatus 25 
for the finite element method according to the present 
Invention will be explained in detail. 

Referring to Rg. 17, the data forming apparatus for 
the finite element method according to this embodiment 
Includes computer 100, fundamental shape storage 30 
device 200, model data storage device 400, analysis 
result storage device 500, display device 600, input 
device 700, printer 800 and recording medium 900. 
Recording medium 900 is a computer readable record- 
ing medium such as a CD-ROM, a magnetic disk, a 35 
semiconductor memory or the like. Programs recorded 
in recording medium 900 are read out by computer 100 
to control the operation of computer 100 and cause 
computer 100 to operate as part data forming means 
101, module data forming means 102. analysis condi- 40 
tion data forming means 103 and data conversion 
means 104. Further, these programs also cause com- 
puter 100 to operate as finite element method analyzer 
105. 

Some fundamental shape data 201 are registered 4S 
in fundamental shape storage device 200 such as a 
magnetic disk or the like in advance. Each of fundamen- 
tal shape data 201 represents a rectangular parallelopi- 
ped, a triangular prism, a cylindrical (column), a 
trapezoidal hexahedron, a notched pillar, a square cylin- so 
der or the like as shown in Figs. 4 through 10. 

According to an indication by a user, part data form- 
ing means 101 appends the actual dimension, the phys- 
ical constants (for example, Young's modulus, Poisson's 
ratio, coefficient of linear expansion, and the like) and 55 
the mesh divisional number to fundamental shapes reg- 
istered in fundamental shape storage device 200, 
thereby forming part data corresponding to the funda- 



mental shapes. The processing of part data forming 
means 101 is shown in Rg. 18. When there Is a selec- 
tion input of a fundamental shape from input device 700, 
part data forming means 1 01 reads out the correspond- 
ing fundamental shape data 201 from fundamental 
shape storage device 200 to display an image of the 
fundamental shape on a screen of display device 600 
(S11). Subsequently, part data forming means 101 
receives the input of the dimension and the physical 
constants corresponding to the fundamental shape from 
input device 700 (812), and further, the input of mesh 
division information as to how the fundamental shape is 
mesh-divided (SI 3). The items of Indicating the dimen- 
sion and the divisional number for each fundamental 
shape may be the same as shown in Fig. 1 1 . Part data 
forming means 101 makes part data obtained, for 
example, by collecting the part name, the selected fun- 
damental shape and the received dimension, physical 
constants and mesh divisional number In a predeter- 
mined format, and stores these part data as one part 
data 301 1 in model data storage device 300 (S14). The 
above operation is repeated for every fundamental 
shapes which are selected by the user (S15). Accord- 
ingly, part data group 301 which consists of required 
number of part data 3011 is generated in model data 
storage device 300. 

Next, module data forming means 102 forms a data 
of module which consist of combination of pairs of parts 
according to a user's input in which the user selects the 
pairs and inputs relative distance between the parts in 
each of the pairs. Fig. 20 is a flowchart showing the 
operation of module data forming means 102. 

When there is a selection Input of a part serving as 
a starting point for the assembly from input device 700, 
module data forming means 102 reads out the corre- 
sponding part data from part data group 301. and dis- 
plays an image of the shape of the part on the screen of 
display device 600 (S21). Subsequently, module data 
forming means 102 receives a selection input of a part 
to be assembled from Input device 700, and reads out 
the corresponding part data from part data group 301 to 
display an image of the shape of the part on the screen 
of display device 600 (S22). Subsequentiy. module data 
forming means 102 displays the screen as shown In Rg. 
14, receives a part having the apex serving as a base 
point, the apex, an apex serving as an origin of a part to 
be assembled, the rotational angle of the latter part, the 
relative distance between both the apexes, and the like 
from Input device 700 (S23), and displays an image of 
an assembly status on the screen of display device 600 
as shown in Fig. 13 (824). If there is input of acceptance 
from input device 700 (YES: S25), module data forming 
means 102 generates module data by grouping In a pre- 
determined format a module name, the name of a part 
to be assembled, tiie apex serving as the origin, rota- 
tional angle of the part, the name of a part serving as a 
base point for the assembly, and the relative distance 
between the apex serving as the base point and the ori- 
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gin, and stores the module data as one module data 
3012 in model data storage device 300 (S26). If there is 
cancel input (not Input of acceptance] (NO: S25), the 
processing returns to step S23 to receive the relative 
distance, and the like again. 

The above operation Is repeated for every assem- 
bly parts selected by the user (S27), whereby an image 
of the final shape Is displayed on display device 600, 
and a module data group 302 which consists of the 
required number of nrKXiule data 3021 is formed. 

Next, analysis condition data forming means 103 
forms analysis condition data 303 according to an anal- 
ysis condition which is indicated by the user, and stores 
the analysis condition data 303 into model data storage 
device 300. Fig. 22 shows an example of the analysis 
condition data 303. A temperature load condition con- 
tains an initial temperature and a reference tempera- 
ture. A constraint contains a part name to be 
constrained and information on a point constraint, a line 
constraint and a surface constraint. For example, when 
a surface is constrained, a surface number, a constraint 
direction, and the like are set. A load condition contains 
a part name to be constrained and information on a 
point load, a line load and a surface load. For example, 
when a load is applied to a point, the load condition con- 
tains a point number and load values for x-. y-. and z- 
directions. Further, when a load is applied to a surfece, 
the load condition contains a surface number and toad 
values for x-, y- and z- directions. 

When model data 304 comprising part data group 
301 , module data group 302 and analysis condition data 
303 are completed in nruxJel data storage device 300 as 
explained above, data conversion means 104 converts 
model data 304 to analysis model 401 for finite element 
method analyzer 1 05. Rg. 23 shows an exanrple of the 
processing thereof. 

First, data conversion means 104 read out model 
data 304 from model data storage device 300 (S31). 
Subsequentiy, data conversion means 1 04 divides the 
fundamental shape having the dimension by tiie mesh 
divisional number for every individual part data 301 1 to 
create meshes, and generates the data consisting of 
the elements and nodal points of the part (S32). The 
positions of each nodal point, and the like are given in a 
coordinate system containing a prescribed apex of the 
part as its origin. Subsequently, data conversion means 
104 converts the data of the elements and nodal points 
of all the parts to data in one coordinate system on the 
basis of tiie relative position between the parts indicated 
by module data group 302 (S33). Subsequently, data 
conversion means 104 converts tine constraint condition 
in analysis condition data 303 to constraint condition 
data for every nodal points (S34), and also the toad con- 
dition in analysis condition data 303 to load condition for 
every nodal points (S35). Subsequentiy, data conver- 
sion means 104 generates constraint condition data for 
connecting nodal points which are non-connected 
between neighboring parts (S36). Thereafter, data con- 



version means 104 outputs analysis model 401 com- 
prising nodal point data, element data, material data, 
constraint data, temperature data, load data and con- 
straint condition data as shown in Fig. 24 to analysis 
5 model storage device 400 (S37). 

The generation of the constraint condition data at 
the step S36 will be explained in more detail. 

Tlie constraint equation on the displacement of a 
structure is represented by the following equatbn (2): 

10 

^/=Z«ctft^q/ (2) 

15 

Uj ; dependent nodal point displacement 
Uqj ; independent nodal point displacement 
T\ ; number of independent nodal points 
Oqj ; coefficient of independent nodal point dis- 
20 placements 

TTie following equation (3) Is used to apply the at>ove 
equation (2) to a connecting condition of nodal points 
between isoparametric three-dimensional 8 nodal point 
25 elements: 

8 

Xi^Y^NqjXqj (3), 

8 

35 8 

/-I 



40 A/<7l«(1+4)(1+tl)(1+C)/8 
/VQ2=(1-4)(1+Ti)(1+0/8 
Nq3=(1-y(1-Ti)(1+C)/8 

45 

A/<74=(1H)(1-tl)(1+C)/8 
A/Q5«(1H)(1+n)(1-C)/8 
50 A/Q6=(1-y(1+Ti)(1-C)/8 
A/(77=(1-^)(1-Ti)(1-0/8 
A/Q8=:(1+«(1-Ti)(1-C)/8 

55 

where 

Xj. Yj, Zf : dependent nodal point coordinate 
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Xqj, Yqj, Zqj : independent nodal point coordinate 
Nqj : three-dimensional one-order shape function of 
Lagrange. 

Such relation ship as shown in fig. 25 is now s 
assumed as an element connection relationship 
between neighboring parts as in the case shown in Fig. 
3. Here, q1 to q8 represent independent nodal points, 
and i represents a dependent nodal point, t r\, C of 
equation (3) which satisfy the element arrangement io 
shown in Fig. 25 are determined (Fig. 26). 4, ti, ^ thus 
calculated are substituted for equation (3) to calculate 
Nqj 0=1 to 8). Nqj are treated as Oqj of equation (2) (j=1 
to 8). and set as the constraint condition data represent- 
ing the element connection. Such constraint condition is 
data are detennined for every non-connected nodal 
points. 

Analysis model 401 thus generated is thereafter 
input to finite element method analyzer 105 and proc- 
essed therein, and analysis result 501 Is output to anal- zo 
ysis result storage device 500 or printer 800. In finite 
element method analyzer 105, the calculation of the 
nodal point displacement or the nodal point temperature 
when non-connected nodal points are contained is per- 
formed using approximating calculation means 106. zs 
Approximating calculation means 106 approximately 
calculates the non-connected nodal point displacement 
or the non-connected nodal point temperature from the 
nodal point displacement or the nodal point temperature 
of the neighboring structural elements on the basis of 30 
the constraint equation. 

As described above, according to the present 
invention, the following effects can be achieved. 

According to the data forming method for the finite 
element method of the present invention, the work of 35 
forming the analysis target shape and for the mesh 
dividing can be performed in a short time. This is 
because it is sufficient for a user to pile up structural ele- 
ments irrespective coincidence or non-coincidence of 
nodal points after inputting a dimensional value for each 40 
of some structural elements and performing mesh gen- 
eration. 

According to the calculation method for the finite 
element method of the present invention, the solution 
can be obtained even when nodal points are not coinci- 45 
dent between elements. This is because when struc- 
tural data that is analysis target in the finite element 
method contain non-connected nodal points, the non- 
connected nodal point displacement or the non-con- 
nected nodal point temperature can be approximated so 
from the nodal point displacement or the nodal point 
temperature of the neighboring structural elements by 
the constraint equation. 

According to the data forming apparatus for the 
finite element method of the present invention, the anal- ss 
ysis model can be formed in a short time. This is 
because after a user appends actual dimension, physi- 
cal constants and mesh divisional number for some fun- 



damental shapes to form parts, inputs the relative 
positions between parts to assemble plural parts on the 
screen of the display device, and inputs analysis condi- 
tions such as the constraint condition, the load condi- 
tion, and the like, the analysis model containing the 
constraint equation for connecting nodal points which 
are non-connected between neighboring parts can be 
automatically generated. 

Claims 

1 . A method for forming structural data to be analyzed 
by a finite element method which comprises steps 
of: 

providing structural elements having various 
shapes; 

inputting dimension for each of said structural 
elements; 

applying mesh generation to each of said struc- 
tural elements; and 

piling up sakJ structural elements in-espective 
of coincidence or non-coincidence of nodal 
points. 

2. A calculation method based on finite element 
method for calculating nodal point displacement or 
nodal point temperature in case that structural data 
to be analyzed by finite element method contain 
non-connected nodal points, which comprises a 
step of calculating the displacements or tempera- 
tures of non-connected nodal points of a structural 
element by using the displacements or tempera- 
tures of nodal points of the neighboring structural 
elements on the basis of constraint equations. 

3. The method as set forth In claim 2. which further 
comprises a step of generating said constraint 
equations by using an isoparametric shape func- 
tion. 

4. An apparatus for generating input data to be ana- 
lyzed by a finite element metiiod, which comprises: 

a storage device which stores fundamental 
shapes; 

part data forming means for setting the actual 
dimension, the physical constants and the 
mesh dividing number for said fundamental 
shapes in accordance with a user's instruction, 
thereby forming the part data for the fundamen- 
tal shapes; 

module data forming means for forming the 
data of a module on the basis of a pair of parts 
indicated by tiie user and the relative position 
between the parts in said pair, said module 
being constructed by a combination of the 
paired parts; 
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analysis condition data forming means for between said parts, 

forming analysis condition data according to an 
analysis condition indicated by the user; and 
data converting means for inputting model data 
comprising said part data, said module data s 
and said analysis condition data, and forming 
an analysis model conprising structural data, 
converted analysis condition data and con- 
straint equations, wherein said structural data 
comprises elements data and nodal points io 
data which are obtained by building up said 
parts in accordance with sard relative positions 
and then dividing the shapes of said parts by 
said mesh division number, each components 
of said converted analysis condition data corre- is 
spends to each nodal points, and said con- 
straint equations connect non-connected nodal 
points between said parts. 

The apparatus as set forth in claim 4. which further 20 
comprises means for generating said constraint 
equations by using an isoparametric shape func- 
tion. 

A computer program product comprising a compu- ss 
ter useable medium having computer program logic 
stored therein, said computer program logic com- 
prising: 

part data forming means for setting the actual 30 
dimension, the physical constants and the 
mesh dividing number for said fundamental 
shapes stored in a storage device in accord- 
ance with a user's instruction, thereby forming 
the part data for the fundamental shapes: 35 
module data forming means for forming the 
data of a module on the basis of a pair of parts 
Indicated by the user and the relative position 
between the parts in said pair, said module 
being constructed by a combination of the 40 
paired parts; 

analysis condition data forming means for 
forming analysis condition data according to an 
analysis condition indicated by the user; and 
data converting means for inputting model data 45 
comprising said part data, said module data 
and said analysis condition data, and forming 
an analysis model conprising structural data, 
converted analysis condition data and con- 
straint equations, wherein said structural data so 
comprises elements data and nodal points 
data which are obtained by building up said 
parts in accordance with said relative positions 
and then dividing the shapes of said parts by 
said mesh division number, each components ss 
of said converted analysis condition data corre- 
sponds to each nodal points, said constraint 
equations connect non-connected nodal points 
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